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We study thermodynamics of SU(3) gauge theory at fixed scales on the lattice, where we vary 
temperature by changing the temporal lattice size A', = (Ta t ) . In the fixed scale approach, finite 
temperature simulations are performed on common lattice spacings and spatial volumes. Con- 
sequently, we can isolate thermal effects in observables from other uncertainties, such as lattice 
artifact, renormalization factor, and spatial volume effect. Furthermore, in the EOS calculations, 
the fixed scale approach is able to reduce computational costs for zero temperature subtraction 
and parameter search to find lines of constant physics, which are demanding in full QCD simu- 
lations. As a test of the approach, we study the thermodynamics of the SU(3) gauge theory on 
isotropic and anisotropic lattices. In addition to the equation of state, we calculate the critical 
temperature and the static quark free energy at a fixed scale. 
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1. Introduction 

Since finite temperature (T > 0) lattice QCD is performed on lattice with a temporal extent 
N t = 1/aT, qualitative calculations at high T may require lower simulation cost than that at T = 0. 
However quantitative systematic studies at T > need huge simulation cost often more than that at 
T = 0. Because such study requires T = calculations at wide range of lattice scale. It is a reason 
why recent large scale thermodynamics calculations are often performed with the Staggered type 
quark formulations [jjj], which needs lower computational cost than that with Wilson type formu- 
lations [||], the domain wall and overlap quarks are all the more costly. To make matters worse, 
some Wilson type quarks sometimes cause some problems at coarse lattice, e.g. nonperturbative 
clover coefficient c$w is reliably determined at a ~ 0. 1 fm or finer [Q], and the domain wall quarks 
encounter with strong residual quark mass effects at coarse lattice p]|. In spite of the difficulties, 
results at T > with the Wilson type quark formulations are desired, because the Staggered type 
quarks suffer from problems of the flavor symmetry violation and the rooted Dirac operators. 

Therefore we propose an alternative fixed scale approach to study thermodynamics of QCD, 
where we vary T by varying the temporal lattice size N t = (Ta t )~ l instead of the conventional fixed 
N t approach. In the fixed scale approach, T = results are common for each N t (T) simulations. It 
may be able to reduce total simulation cost drastically. Furthermore, common parameters (except 
for N t ) enable us to investigate pure thermal effects of observables without obstacles coming from 
changing lattice spacing and spatial volume effects. 

In this report, we test the approach in the SU(3) gauge theory on isotropic and anisotropic lat- 



tices. Our lattice action and some details of the EOS calculation are given in Sect J2.1| and Sect.[2.2 



Results of EOS are presented in Sects. 2.3 and 2.4. The T c and the static quark free energy are 



discussed in Sect.|| and^. We conclude in the last section. 

2. Equation of state 

2.1 T-integration method 

In the fixed scale approach, to calculate the pressure non-perturbatively, we propose a new 
method, "the T-integral method" |Q] : 

based on another thermodynamic relation valid at vanishing chemical potential: 

T ±(P_\e-3p 
dT \T 4 J T 4 

The initial temperature Tq is chosen such that p(Tq) rj 0. Calculation of e — 3p requires the beta 
functions just at the simulation point, but no further Karsch coefficients are necessary. Since T is 
restricted to have discrete values, we need to make an interpolation of (e — 3p) /T 4 with respect to 
T. 

Since the coupling parameters are common to all temperatures, our fixed scale approach with 
the T-integral method has several advantages over the conventional approach; (i) T = subtractions 
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can be done by a common T = simulation, (ii) the condition to follow the LCP is obviously 
satisfied, and (iii) the lattice scale as well as beta functions are required only at the simulation point. 
As a result of these, the computational cost needed for T = simulations is reduced drastically. 

When we adopt coupling parameters from T = spectrum studies, the values of N t around 
and below the critical temperature T c are much larger than those used in conventional fixed N t 
studies. For example, at a m 0.07 fm, T ~ 175 MeV is achieved by N t ~ 16. Therefore, for 
thermodynamic quantities around and below T c , we can largely reduce the lattice artifacts over the 
conventional approach, with much smaller total computational cost. This is also a good news for 
phenomenological applications of the EOS, since the temperature achieved in the relativistic heavy 
ion collision at RHIC and LHC will be at most up to a few times the T c 0]. We note here that, as T 
increases, N t becomes small and hence the lattice artifact increases. Therefore, our approach is not 
suitable for studying how the EOS approaches the Stephan-Boltzmann value in the high T limit. 



2.2 Lattice action 

We study the SU(3) gauge theory with the standard plaquette gauge action on an anisotropic 
lattice with the spatial (temporal) lattice size and scale N s (N t ) and a s (a t ), respectively. The lattice 
action is given by 



3 

r 

X i=l 



1 - -ReTrIfo(jc) 



X i>j=l 



l--ReTrC/y(jc) 



def. 



(2.3) 
(2.4) 



where U jlv (x) is the plaquette in the /iv plane and j3 and £0 are the bare lattice gauge coupling and 
bare anisotropy parameters. The trace anomaly is calculated by 



e-3p 




1 



(2.5) 
(2.6) 



with £ = a s /a t the renormalized anisotropy. a s (dfi/da s ) is the beta function, d^o/df} vanishes on 
isotropic lattices. 



2.3 EOS on isotropic lattice 

Our simulation parameters are listed in Table |]. On isotropic lattices, we calculate EOS on 
three lattices to study the volume and lattice spacing dependences. The ranges of N t correspond 
to T = 210-700 MeV for the sets il and i2, and T = 220-730 MeV for i3, to be compared with 
T c ~ 290 MeV. The set a2 will be discussed later. The T = subtraction is performed with N t = 16 
for il and i2, and with N t = 22 for i3. We generate up to a few millions configurations using the 
pseudo-heat-bath algorithm. Statistical errors are estimated by the Jackknife analysis, appropriate 
bin sizes, which strongly depend on T. Typically, bin size of a few thousands configurations are 
necessary near T c , while a few hundreds are sufficient off the transition region. 
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set 


P 




A(s 




r /a s 


a s [fm] 


L[fm] 


<a(<i,g 2 /da) 


il 


6.0 


1 


16 


3-10 


5.35(+2) 


0.093 


1.5 


-0.098172 


i2 


6.0 


1 


24 


3-10 


5.35(+2) 


0.093 


2.2 


-0.098172 


i3 


6.2 


1 


22 


4-13 


7.37(3) 


0.068 


1.5 


-0.112127 


a2 


6.1 


4 


20 


8-34 


5.140(32) 


0.097 


1.9 


-0.10704 



Table 1: Simulation parameters on isotropic and anisotropic lattices. On isotropic lattices, we adopt ro/aof 
[||], and the beta function of [||]. Anisotropic ro/a s is from [[To|], and the beta function is calculated in [Q], 
Lattice scale a 4 and lattice size L = N s a s are calculated with rg = 0.5 fm. 
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Figure 1: (Left) Trace anomaly on isotropic lattices. The dotted lines are natural cubic spline interpolations. 
Horizontal errors due to the lattice scale are smaller than the symbols. (Right) The energy density and the 
pressure on isotropic lattices. 



Figure [I] (Left) shows (e — 3p)/T 4 . Dotted lines in the figure are the natural cubic spline 
interpolations. At and below T c , lattice size dependence is visible between the sets il (L ss 1.5 fm) 
and i2 (2.2 fm). On the other hand, the lattice spacing dependence is small between il {a 0.093 
fm) and i3 (0.068 fm). At higher T, (e — 3p)/T 4 on three lattices show good agreement. The 
integration of (2.1) is performed numerically using the natural spline interpolations shown in Fig.|l] 
(Left). For the initial temperature To of the integration, we linearly extrapolate the (e — 3p) /T A data 
at a few lowest T's because the values of (e — 3p)/T at our lowest T are not exactly zero. In this 
study, we commonly take To = 150 MeV as the initial temperature which satisfies (e — 3p) /T 4 = 0, 
and estimate the integration from To to the lowest T by the area of the triangle. Statistical errors 
for the results of integration are estimated by a Jackknife analysis [11]. Note that the error in the 
lattice scale do not affect the dimension less quantity p/T 4 . 

In Fig.|l] (Right), we summarize the results of EOS on isotropic lattices. e/T 4 is calculated 
combining the results of p/T 4 and (e — 3p)/T 4 . Since the lattice parameter dependence of (e — 
3p) /T 4 is small except for the vicinity of T c , we find that EOS has a similar shape except for the 
vicinity of T c . Near and below T c , we observe a sizable finite volume effect between L sa 1.5 fm 
and 2.2 fm, while the lattice spacing effects are not so. At large T, we note a slight tendency that 
p and £ decrease as the lattice size becomes larger and the lattice spacing becomes smaller. Our 
results are qualitatively consistent with the previous results by the conventional fixed N t method 
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Figure 2: (Left) The trace anomaly on the anisotropic lattice (a2) and the isotropic lattice (i2). Dotted lines 
show the cubic spline interpolation. (Right) The EOS on the anisotropic lattice. 

but with much reduced lattice artifacts around T c due to much larger there. 

2.4 EOS on anisotropic lattice 

The anisotropic lattice with the temporal lattice finer than the spatial one is expected to improve 
the resolution of T without much increasing the computational cost. To further test the systematic 
error due to the resolution of T, we perform the study with the T -integral method on an anisotropic 
lattice with the renormalized anisotropy ^ = 4. The simulation parameters are given as the set a2 
in Table [j], which are the same as those adopted in JTo|]. We vary N t = 34-8 corresponding to 
T = 240-1010 MeV. The T = subtraction is performed with N t = 80. We generate up to a few 
millions configurations. The beta function a s (dp /da s )p with = 4 has calculated in our paper [j|], 
and its value at our simulation point is given in Table [jj. For (d^o/df})? we adopt the result of [p"2|]. 

In Fig. |2| (Left), we compare the trace anomaly obtained on the anisotropic lattice with that on 
the isotropic lattice with similar a s and L (the set i2). We find that the results are generally consistent 
with each other except for around T c . We note a systematic tendency that the trace anomaly on the 
anisotropic lattice is slightly lower than that on the isotropic lattice. According to this tendency, 
the pressure on the anisotropic lattice is slightly smaller than that on the isotropic lattice at high T. 
The tendency may be understood by the smaller lattice artifact due to the temporal lattice spacing 
on anisotropic lattices, since lattice artifacts due to temporal lattice spacing are larger than that by 
the spatial lattice spacings in thermodynamic quantities [|T3|]. Finally, we summarize our results 
of EOS on the anisotropic lattice in Fig. ^| (Right). We find that they are consistent with those on 
isotropic lattices. 

3. Transition temperature 

Here we consider a possibility to compute the T c in the fixed scale approach. T c is determined 
by studying temperature dependence of order parameters. Strictly speaking, such temperature de- 
pendence should be separated from other effects, such as renormalization, lattice artifacts, and 
spatial volume dependence. In the fixed scale approach, we can easily isolate the thermal effect on 
the observables. On the other hand, the resolution of T is restricted by descrete N t . 
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set 


P 






N t 


ro/a s 


a s [fm] 


L[fm] 


a2-l 


6.1 


4 


20 


26-30 


5.140(32) 


0.097 


1.9 


a2-2 


6.1 


4 


30 


26-30 


5.140(32) 


0.097 


2.9 


a2-3 


6.1 


4 


40 


26-30 


5.140(32) 


0.097 


3.9 



Table 2: Simulation parameters on anisotropic lattices to study volume dependence of the susceptibility of 
the Polyakov loop. 




Figure 3: (Left, Center) Polyakov loop and its susceptibility on the anisotropic lattice. (Right) Static quark 
free energy on the anisotropic lattice. 

We calculate the Polyakov loop and its susceptibility on the anisotropic lattice. In addition 
to the lattices prepared for the EOS calculation, we generated different spatial volume lattices to 
study its finite size scaling. The parameters we adopted are listed in Table ||[ Since the SU(3) gauge 
theory has the global Z(3) symmetry, we calculate the real part of Z(3)-rotated Polyakov loop and 
its susceptibility. Our results are shown in the left and center pannels of Fig.||. Unlike the case of 
conventional fixed N t studies, the renormalization factor is common to all temperatures. From the 
susceptibility data we can find that the transition point locates at N t 28 corresponds to T c = 280- 
300 MeV from the vq scale setting. We also find that the peak height of the susceptibility increases 
with increasing the system volume, in accordance with the 1st order nature of the transition. 

4. Static quark free energy 

Finally, we study the static quark free energy V(r) in the fixed scale approach. In conventional 
fixed N t studies, the additive renormalization constant for V (r) is diffrent for each T because the 
lattice spacing is different. Assuming that the short distance physics is independent of T, the 
constant term is conventionally adjusted by hand such that V(r) around the smallest r coinside 
with each other. 

In the fixed scale approch, on the other hand, the common lattice spacing for all T implies that 
the constant term in V(r) should be common too. Therefore, we can purely study the T effects 
without adjusting the constant term. In Fig.|| (Right) we show our results of V(r) for the color 
singlet channel on the anisotropic lattice, without adjusting the constant term. The solid curve in 
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the figure is a fit result of T = potential. We find that V(r) at different temperatures converge to 
a common curve at short distances. We thus have confirmed the expectation that the short distance 
physics is independent of T. 

5. Conclusions 

We proposed a fixed scale approach to study the QCD thermodynamics on the lattice. In 
this approach, T is varied by changing the temporal lattice size N t at a fixed lattice scale. To test 
the method, we applied it to the SU(3) gauge theory on isotropic and anisotropic lattices. We 
found that the T-integral method to calculate the EOS works quite well. The main advantage 
of our approach is that the computational cost for T = simulations, which are the most time 
consuming calculations in the conventional fixed N t approaches, can be drastically reduced. We 
may even borrow configurations of existing high precision simulations at T = 0. The approach 
is applicable to QCD with dynamical quarks too. We are currently investigating EOS in 2 + 1 
flavor QCD with non-perturbatively improved Wilson quarks, using the configurations by the CP- 
PACS/JLQCD Collaboration [[14]]. With these fine lattices, the lattice artifacts around T c are much 
smaller than the conventional fixed Nt approaches. We are further planning to use the PACS-CS 
configurations just at the physical point |fj]. 

TU thanks H. Matusufuru for helpful discussions and comments. The simulations have been 
performed on supercomputers at RCNP, Osaka University and YITP, Kyoto University. This work 
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